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SUR  LES  DÉFORMATIONS  ET  L’EXTINCTION  DES  ONDES  AÉRIENNES,  ISOLÉES 
OU  PÉRIODIQUES,  PROPAGÉES  A L’INTÉRIEUR  DE  TUYAUX  DE  CONDUITE 
SANS  EAU,  DE  LONGUEUR  INDÉFINIE; 

Par  M.  J.  BOU£SINESQ. 

I.  La  propagation  du  son,  le  long  des  tuyaux  de  conduite  pleins 
d’air,  a été  étudiée  théoriquement,  dans  l’hypothèse  de  mouve- 
ments pendulaires  d’assez  faible  amplitude  pour  permettre  la  ré- 
duction des  équations  à la  forme  linéaire,  d’abord  par  M.  von 
Helmholtz  (1),  en  négligeant  les  échanges  alternatifs  de  chaleur 
entre  la  paroi  et  la  couche  gazeuse  contiguë,  puis  par  Kirchhoff  (2) 
en  tenant  compte  de  ces  échanges.  Les  formules  auxquelles  ils 
sont  parvenus,  différentes  seulement  par  la  signification  du  coef- 
ficient au  moyen  duquel  s’y  exprime  l’influence  extinctrice  et 
retardatrice  de  la  paroi,  ont  été  récemment  employées  par 
MM.  Violle  et  Vauthier,  dans  un  Mémoire  des  Annales  de 
Chimie  et  de  Physique  (mars  1890,  6e  série,  t.  XIX)  où  ceux-ci 
ont  publié  et  discuté  leurs  remarquables  observations,  faites  à 
Grenoble,  en  1 885,  sur  un  tuyau  de  conduite  neuf  de  om,  70  de 
diamètre  et  près  de  i3km  de  longueur,  parcouru  plusieurs ‘fois  par 
chaque  onde  que  réfléchissaient  successivement,  aux  deux  extré- 
mités, des  parois  planes  y fermant  le  tuyau  perpendiculairement 
à l’axe.  Mais  ces  observations  n’ont  guère  été  relatives  qu’à  des 
ébranlements  non  périodiques,  d’une  durée  totale  restreinte  au- 
près de  chaque  section,  et  dont  les  plus  caractéristiques  consi- 
staient en  ondes  condensées  produites,  comme  dans  de  mémo- 
rables expériences  antérieures  de  Régnault,  par  des  coups  de 
pistolet  tirés  à une  extrémité  du  tuyau,  devant  un  orifice  percé 
au  centre  de  son  fond  plat  et  que  l’on  bouchait  aussitôt  après.  Or, 
malgré  quelques  systèmes  de  mouvements  vibratoires  qui,  surtout 
dans  les  premiers  instants  de  la  propagation,  se  superposent  à 
fonde  générale,  ou  la  sillonnent  transversalement,  et  lui  per- 
mettent d’affecter  l’organe  de  l’ouïe,  de  telles  intumescences 

(')  Verhandl.  d.  natur.  Vereins  zu  Heidelberg , t.  III,  p.  16;  i863. 

(a)  Poggendorff’s  Annalen,  l.  CXXXIV,  p.  177;  1868. 
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où  l’air  est  condensé  presque  d’un  bout  à l’autre,  diffèrent  pro- 
fondément de  la  suite  périodique  d’ondes,  à égales  dilatations 
et  condensations  alternatives,  supposée  par  les  formules  de 
MM.  HelmhoJtz  et  Kirchhoff.  Aussi  MM.  Violle  et  Vauthier  ont-ils 
trouvé  que  celles-ci  indiquaient,  pour  un  parcours  donné,  un 
affaiblissement  de  la  pression  incomparablement  plus  rapide  que 
le  vrai  et,  au  contraire,  un  ralentissement  de  la  propagation 
presque  négligeable  vis-à-vis  des  petits  ralentissements  réels,  ob- 
servés déjà  par  Régnault. 

La  principale  raison  d’aussi  graves  désaccords  ne  se  trouve  sans 
doute  pas  dans  l’assimilation,  qui  doit  être  permise  approximati- 
vement ou  comme  premier  aperçu,  d’une  intumescence  isolée, 
d’une  certaine  longueur  apparente,  à la  moitié  convexe  d’une  onde 
d’un  système  périodique  où  chaque  demi-ondulation  aurait  cette 
longueur.  Mais  elle  est  évidemment  dans  la  confusion  établie  entre 
l’intumescence  entière  cheminant  le  long  du  tuyau  et  les  ondes 
sonores  successives,  beaucoup  plus  courtes,  qui  la  sillonnent, 
auxquelles  seules  conviennent  le  ralentissement  et  le  coefficient 
d’extinction  calculés.  Ceux-ci,  en  effet,  d’après  les  formules  théo- 
riques, sont,  l’un,  c’est-à-dire  le  ralentissement,  en  raison  directe 
et,  l’autre,  le  coefficient  d’extinction,  en  raison  inverse,  de  la  racine 
carrée  de  la  durée  de  vibration,  ou  de  la  racine  carrée  de  la  lon- 
gueur d’onde;  et  ils  deviendraient,  par  suite,  le  premier,  bien  plus 
sensible,  mais  le  second,  bien  plus  faible,  si  on  les  évaluait  pour 
des  ondes  aussi  longues  que  l’intumescence  entière,  à laquelle 
MM.  Violle  et  Vauthier  ont  effectivement  reconnu  beaucoup  plus 
de  longévité  ou  de  persistance  qu’aux  sons  accessoires  qui  l’accom- 
pagnent au  début. 

II.  Ainsi,  les  désaccords  sont  plus  apparents  que  réels,  malgré 
Ja  portée  restreinte  ou  le  défaut  de  généralité  des  formules.  Néan- 
moins, il  y a évidemment  lieu  de  reprendre  le  problème  au  point 
de  vue  théorique,  en  attribuant,  aux  condensations  et  aux  vitesses 
successivement  produites  sur  une  même  section  quelconque,  des 
expressions  non  plus  pendulaires  par  rapport  au  temps,  mais 
affectées  d’une  fonction  arbitraire,  propre  à représenter  les  diverses 
successions  possibles  de  pressions^ou  de  déplacements  directement 
réalisables  près  de  l’entrée  du  tuyau.  Tel  est  l’objet  du  présent 
travail. 
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Une  question  assez  analogue  d’ondes  liquides,  mais  plus  simple, 
celle  de  l’extinction  graduelle  d’une  intumescence,  d’une  onde 
solitaire  par  exemple,  le  long  d’un  canal,  par  les  frottements  pro- 
duits surtout  près  du  fond  et  des  bords,  peut  suggérer  la  méthode 
à suivre  : je  l’ai  traitée  au  § II  (nos  9 à 19)  d’un  Complément  à 
ma  Théorie  des  eaux  courantes,  publié  en  octobre  1878  dans  le 
Journal  de  Mathématiques  pures  et  appliquées  (t.  IV,  voir  les 
p.  347  à 366).  Là  aussi,  les  principales  résistances  au  mouvement 
sont  localisées  dans  une  mince  couche  contiguë  aux  parois,  où  les 
vitesses  varient  très  rapidement,  sur  une  épaisseur  insensible, 
depuis  la  valeur  zéro  maintenue  constamment  sur  la  face  exté- 
rieure, jusqu’aux  valeurs  complètes  relatives  à l’intérieur  de  la 
masse  fluide  et  réalisées  déjà  sur  l’autre  face  de  la  couche,  valeurs 
pouvant  être  une  fonction  arbitraire  donnée  du  temps.  Les  vitesses 
et  les  frottements,  dans  l’épaisseur  de  cette  couche,  s’expriment 
par  des  intégrales  définies  simples,  prises  de  zéro  à l’infini,  mais 
dont  malheureusement  la  sommation,  possible  à la  rigueur  après 
spécification  de  la  (onction  arbitraire,  ne  paraît  eflectuable  sous 
forme  finie  que  dans  peu  de  cas,  parmi  lesquels  le  plus  simple  peut- 
être  est  justement  celui  de  mouvements  pendulaires.  D’ailleurs, 
le  serait-elle  toujours,  que  l’utilisation  de  ses  résultats  exigerait 
encore  une  connaissance  assez  précise,  théorique  ou  expérimen- 
tale, de  la  forme  affectée  par  chaque  variété  d’ondes,  forme  dont 
dépendent  les  circonstances  cherchées  de  la  propagation.  De  là, 
une  grande  difficulté  pour  calculer  l’affaiblissement  graduel  des 
intumescences.  Mais,  du  moins,  quand  il  s’agit  d’une  onde  soli- 
taire, dont  le  profil,  éminemment  stable  pour  chacune  des  valeurs 
successives  décroissantes  de  l’énergie  totale,  se  trouve  représentée 
par  une  fonction  toujours  la  même  d’une  variable  à expression 
linéaire  par  rapport  à l’abscisse  et  monôme  par  rapport  à la  hauteur 
ou  à l’énergie  totale  de  l’onde,  le  lent  décroissement  de  cette 
hauteur  ou  de  cette  énergie  se  détermine  facilement,  à l’exception 
d’un  facteur  constant  dans  le  calcul  duquel  reste  incluse  toute  la 
difficulté  de  la  quadrature  ineffectuée. 

Or,  au  contraire,  les  ondes  aériennes  condensées  se  propageant 
le  long  d’un  tuyau  ne  tendent,  avant  leur  extinction,  vers  aucune 
forme  commune.  Aussi,  quoique  une  double  application  de  la  for- 
mule d’ Analyse,  à intégrale  définie  simple,  utilisée  dans  le  cas  des 
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ondes  liquides,  permette  d’y  exprimer,  d’une  part,  les  frottements 
successifs  de  la  paroi,  avec  les  vitesses  réalisées  dans  la  couche 
conliguë,  d’autre  part,  les  densités  propres  à cette  couche,  et  en- 
suite de  former  les  équations  de  mouvement  des  tranches  gazeuses 
normales  à l’axe,  d’évaluer  les  vitesses  de  propagation  des  divers 
éléments  de  masse  de  l’intumescence,  etc.,  il  n’y  a guère  que  les 
circonstances  générales  de  sa  déformation  et  de  son  extinction  qui 
puissent  facilement  se  reconnaître  sur  les  formules.  Les  princi- 
pales de  ces  circonstances,  bien  d’accord  avec  les  résultats  de  l’ob- 
servation, consistent  dans  un  aplatissement  et  un  allongement 
indéfinis  de  l’onde,  du  moins  après  une  première  période  assez 
courte  où  son  sommet  s’est  rapidement,  s’il  y avait  lieu,  rapproché 
de  son  front,  c’est-à-dire  de  sa  face  antérieure,  qui  reste  d’ailleurs 
très  sensiblement  plane  à toutes  les  distances  de  l’extrémité  de  dé- 
part. 

On  verra  que  la  théorie  indique  de  plus,  entre  le  coefficient 
d’extinction  et  le  ralentissement  causé  par  la  paroi,  une  relation 
approximative,  à laquelle  les  expériences  tant  de  Régnault  que  de 
MM.  Violle  et  Vauthier  apportent  une  certaine  confirmation. 

III.  Avant  d’étudier  l’influence  de  la  perméabilité  calorifique  et 
du  frottement  de  la  paroi  du  tuyau  sur  les  ondes,  censées  venues 
de  loin,  que  propage  la  masse  gazeuse  remplissant  celui-ci,  sup- 
posons nulle,  pour  un  instant,  cette  double  influence,  afin  de  pou- 
voir admettre  partout  la  forme  plane  des  ondes  avec  le  parallélisme 
des  mouvements  à l’axe  (rectiligne),  c’est-à-dire  tout  à la  fois,  la 
parité  absolue,  à chaque  instant,  de  la  température,  de  la  pression 
et  de  la  vitesse,  sur  toute  l’étendue  d’une  section  normale  quel- 
conque du  tuyau,  comme  si  le  milieu  gazeux  était  latéralement  indé- 
fini. Alors,  dans  les  équations  indéfinies  de  mouvement  du  fluide, 

i / «?N.r  dTz  dTy\_  , 

p \ dx  dy  + dz  J ~ U ’ 
i /dT~  «?Nr  , 

p \ dx  "+*  dy  dz  ) V ’ 

I fdTy  dTx  dNz\  , 

p \ dx  dy  dz  ) W ’ 
dp  d.pu  d.pv  d.ow 

dt  dx  dy  dz 
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où  les  composantes  N,  T de  pression  ont  les  formules  respectives, 
à deux  coefficients  de  frottement  intérieur  ou  d’imparfaite  fluidité  e, 
£,  — 2 S, 


(2) 


/ du  dv  dw 
\ dx  ’ dy  ’ dz 


les  composantes  transversales  c,  w de  la  vitesse  seront  nulles,  et 
sa  composante  u , longitudinale,  ou  parallèle  à l’axe  du  tuyau 
choisi  pour  axe  des  x,  dépendra  seulement  de  x et  de  t , ainsi  que 
la  densité  p,  la  température  absolue  T et,  par  suite,  la  pression 
élastique  /?,  proportionnelle  au  produit  oT. 

Appelant  p0 , p0,  T0  la  pression,  la  densité  et  la  température 
constantes  primitives,  posons 


(3) 


P — Po ( i Y ) » T — T0(i-4-t), 


où  y,  t désignent  ainsi  la  condensation  et  réchauffement  relatifs 
dus  au  passage  des  ondes.  On  sait,  par  le  premier  principe  de  la 
Thermodynamique,  que,  si  M désigne  la  masse  d’une  particule 
gazeuse  et  G son  calorique  spécifique  à volume  constant,  la  quan- 
tité de  chaleur  (évaluée,  comme  C,  en  kilogrammètres)  qui 

pénètre  du  dehors  par  sa  surface,  durant  un  instant  dt  où  sa 
densité  et  sa  température  croissent  respectivement  de  dp  et  de  dT , 
a comme  expression,  en  appelant  finalement  m — i le  rapport 


constant  — ^ 
CpT 


(4) 


dQ  = M ( G dT  -h  p d--j  = M ( C dT 


P dp 


= MCT 
= MCT 


r^r 

l.T 


dp 


(±-m±\  = mCTdlos-P-. 


\P 


De  plus,  l’accord  très  satisfaisant  de  la  vitesse  expérimentale 
du  son  avec  sa  formule  théorique  donnée  par  Laplace,  montre 
l’extrême  petitesse  des  échanges  de  chaleur  s’effectuant  dans  les 
sens  parallèles  à la  propagation,  ou  entre  les  couches  gazeuses 
successivement  atteintes  par  les  ondes;  de  sorte  que  l’on  peut 


6 


J.  BOUSSINESQ.  — DÉFORMATION  ET  EXTINCTION 

écrire  ici  6/Q  = o,  et  supposer  la  pression  p proportionnelle  à pm, 
c’est-à-dire  à (i  -4-y)m.  Nous  aurons  donc 

(5)  p — /?0(i  -f-  y)7”  dans  une  onde  latéralement  indéfinie, 


ou  mieux,  en  substituant  à /?0  sa  valeur  en  fonction  de  la  vitesse 
connue  a du  son,  exprimée  par  \/mfo’ 

po«2 


(6) 


P = — (l  H"  Y)' 

m 17 


Enfin,  les  frottements  intérieurs  ou  autres  parties  non  élas- 
tiques des  pressions  seront  assez  faibles,  pour  qu’on  puisse  né- 
gliger, dans  (2),  les  termes  non  linéaires  qui  contiendraient  en 
facteur,  outre  une  dérivée  partielle  des  vitesses  u , v,  w , la  petite 
variation  des  coefficients  e et  e,  — -2e  due  au  changement  relatif  t 
de  la  température,  ou  aussi  à la  condensation  y (si  s,  e,  — 2 s dé- 
pendaient sensiblement  de  la  densité,  ce  qui  ne  paraît  pas  avoir 
lieu),  quantités  évanouissantes  avec  l’amplitude  des  ondes.  En 
d’autres  termes,  les  valeurs  de  s et  de  £1  — 2s  à introduire  seront 
celles  de  ces  coefficients  pour  la  température  et  la  densité  primi- 
tives constantes  p0,  t0.  Et,  de  même,  le  diviseur  p pourra  être  ré- 
duit à p0  dans  les  termes  des  premiers  membres  de  (i)  où  figure- 
ront e,  . 

Cela  posé,  en  observant  que  les  hypothèses  faites  annulent  e, 
w,  v',  w>'  et  rendent  u,  y fonctions  seulement  de  x et  de  les 
valeurs  (2),  (3)  des  N,  T et  de  p,  après  substitution  à p de  son 

expression  (6)  et  à u!  de  ^ -h  u ^ ? réduiront  les  quatre  équa- 
tions (1)  du  mouvement  aux  deux  suivantes  : 


(7) 


du 

~dt  + a 


d y 


du 


r )m~l  -p-  A-  u —j— 
1 dx 


dx 


dt 


du 


d . u y 
dx  dx 


d 2 u 
p0  dx 2 


IV.  A une  première  approximation,  où  sont  négligeables  les 
termes  non  linéaires  en  u et  y,  ainsi  que  celui  où  figure  le  petit  coef- 
ficient s,,  ces  deux  équations  deviennent 

du  dv  dv  du 

Tt+a  Tx=°'  dt+dï=°- 


(8) 
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et  elles  donnent,  par  l’élimination  de  u , 


(9) 


*1 

dt2 


a2S=0- 


L’intégrale  générale  de  cette  dernière  doit  être  spécifiée  pour 
les  ondes  dont  il  s’agit  ici,  propagées  au  sein  d’un  gaz  en  repos. 

Admettant,  par  exemple,  qu’à  une  époque  initiale  très  an- 
cienne, y,  u s’annulassent  pour  toutes  les  valeurs  de  x comprises 
entre  une  abscisse  négative  fort  grande  et  x = cc,  on  peut  con- 
venir de  n’étudier  le  mouvement  ultérieur  qu’en  ces  endroits,  où 
le  fluide  était  primitivement  immobile.  Alors,  si  cp  désigne  une 
fonction,  d’ailleurs  arbitraire,  s’annulant  pour  la  valeur  H-00  de  sa 
variable,  il  vient,  comme  on  sait,  y=zo(x  — at );  après  quoi  la 

seconde  équation  (8),  écrite  — ay  + u)  = o et  intégrée  sous 

la  condition  imposée  que  les  deux  fonctions  y,  u s’annulent  pour  x 
infini  positif,  donne  u — ay,  valeur  rendant  identique  la  première 
(8),  vu  l’expression  de  y.  L’on  a donc,  à une  première  approxi- 
mation, 


(10)  y = cp(a?  — at),  u=-av(x  — at); 

formules  d’où  il  résulte  que  toute  fonction  des  deux  quantités  u , 

/<  » /-»  co 

y,  et  même  toute  intégrale  comme  / y dx  — / cp(a)  dx,  ..., 

d X dx — at 

dépendent  presque  uniquement  de  la  variable  x — at , et  admet- 
tent la  formule  symbolique  approchée  de  différentiation 


(11) 


dt  a dx 


A une  deuxième  et  suffisante  approximation,  les  termes  de  (7) 
qui  ont  pu  être  négligés  dans  (8)  seront  donc  évaluables  par  les 
formules  (10),  (11);  ce  qui  permettra  de  remplacer,  au  besoin, 
dans  ces  termes  ou  dans  leurs  dérivées,  u par  a y et  les  différen- 
tiations relatives  à t par  d’autres  relatives  à x.  Si,  d’ailleurs,  à 
(1  H-  y)w-2  l’on  substitue  non  plus  l’unité,  mais  la  valeur  beaucoup 
plus  exacte  i + (m  — 2)  y,  ces  équations  deviendront 


(12) 


du  „ dy 
dï+a  ~dx 


d_ 
dx  \ 


/ m 


s 1 dy 
■ p—J-a-d- 
po  dx 


dt 


du  d.  y2 

7hc^a~d^ 


= 0. 


Différentiées  la  première,  en  #,  la  deuxième,  en  t , puis  multi- 


8 


J.  B0USS1NESQ.  — DÉFORMATION  ET  EXTINCTION 

pliées  respectivement  par  — i,  4-1  et  ajoutées,  elles  donneront 
comme  équation  en  y,  au  lieu  de  (9),  une  relation  de  la  forme 


(i3) 


= £* 

dt2  dx 2 dx  ’ 


d,f 


d 2 y 


dx  p0  dx 2 


Son  intégration,  pour  les  ondes  dont  il  s’agit,  propagées  vers 
les  x positifs  et  à la  tête  desquelles  s’annule  K,  permet  de  la  rem- 
placer par  l’équation,  du  premier  ordre  en  £, 


04) 


dy  dy  K 

— -h  d — — -H  — O, 

dt  dx  ‘±a 


d’après  une  démonstration  (utile  dans  l’étude  de  plusieurs  sortes 
de  mouvements  ondulatoires),  que  l’on  peut  voir  aux  p.  364* 
et  365*  du  tome  11  de  mon  Cours  d’ Analyse  infinitésimale . En 
substituant  à K son  expression,  il  viendra  donc,  pour  déterminer 
les  déformations  successives  de  l’onde,  l’équation 


05) 


dy  dy  m- 1- 1 d.  y2  £t  d2  y 
dt  dx  4 dx  2p0  dx 2 


V.  Cette  équation  implique  la  conservation,  à toute  époque,  de 
la  somme  f y dx,  prise  sur  toute  la  longueur  de  l’onde  ou  entre  les 
limites  ^ = ztoo,  et  elle  entraînerait  aussi  l’invariabilité  de  l’in- 
tégrale f y 2dx  entre  les  mêmes  limites,  sans  l’imperfection  de  la 
fluidité,  c’est-à-dire  si  l’on  avait  = o.  En  effet,  l’équation  ( 1 5 ) , 
multipliée  soit  par  dx , soit  par  2y dx,  et  intégrée  de  x = — co 

à x — oo:  en  observant,  d’une  part,  que  2 y égale  (f  y3)  et  que 


d2  y , , d2. y2  dy2  ,,  dy 

2Ï  555  eSale  Tfo*  - 2 dk’  d autre  Part’  tIue  Tx 

tête  et  à la  queue  de  l’onde,  donne  simplement 


s’annulent  à la 


(16) 


d_ 

dt 


f*—- 


Les  deux  sommes  / y dx,  f y 2 dx  sont  donc,  la  première,  inva- 
riable, mais,  la  deuxième,  sans  cesse  décroissante,  sauf  quand  e,  =0. 

C’est  ce  qu’on  aurait  pu  prévoir  : car,  <r  désignant  l’aire  de  la  sec- 
tion du  tuyau,  le  produit  de  la  première  par  le  facteur  constant  p0<T, 
savoir  p0(jfydx}  ou  f (p — p0fidx,  représente  évidemment  la 
masse  totale  de  l’intumescence,  quantité  de  matière  gazeuse  ini- 
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tialement  introduite  à l’entrée  du  tuyau,  en  sus  de  l’air  qu’il  con- 
tenait à la  densité  uniforme  p0;  et  l’expression  p0a2<r  f v2dx,  pro- 
portionnelle à fy2dx,  représente  très  sensiblement,  à toute 
époque,  V énergie  totale  de  l’intumescence,  c’est-à-dire  le  travail 
qu’elle  emmagasine,  non  moins  constant  que  sa  masse  quand  on 
néglige  l’imperfection  de  la  fluidité.  Effectivement,  d’une  part, 
V énergie  actuelle  ou  demi-force  vive  de  l’onde,  u2  p(<r  dx), 
c’est-à-dire,  à fort  peu  près,  /^(ay)2  p0ar  dx,  constitue  Ja  moitié  de 
l’expression  considérée  p0a2o-  f y2  dx\  et,  d’autre  part,  son  énergie 
potentielle  ou  de  ressort,  qui  est,  pour  chaque  particule  de  masse  M 
M 

et  de  volume  — > le  travail  que  produirait,  durant  sa  détente,  la 
partie  variable  (seule  active)  p — jO0  de  sa  pression,  savoir 


= (sensiblem.)  — — f my  dy  ou  -Ma2  y2,  é-ïierg.  act.  delà  part., 
Po  d0  2 


sera  de  même,  pour  toute  l’onde,  ^p0 a2  <7/  y 2dx,  ou  l’autre  moitié 
de  l’intégrale  considérée. 


VI.  Imaginons  une  famille  de  plans,  normaux  à l’axe  du  tuyau, 
et  mobiles  de  manière  à intercepter  sans  cesse  entre  eux  une  même 
masse  d’intumescence,  ou  de  manière  que  chacun,  dont#  désignera 
l’abscisse  à l’époque  t,  ait  sans  cesse  devant  lui  une  partie  con- 
stante tu—  / y dx  de  l’intégrale  / y dx,  qui  mesure  propor- 
tionnellement  la  masse  totale  de  l’onde.  Alors  la  vitesse  de 
transport^  de  l’un  quelconque  de  ces  plans  sera  ce  que  l’on 

peut  appeler  la  vitesse  de  propagation  de  l’élément  contigu 
dvs  — — y dx  d’intumescence,  et,  en  la  désignant  par  to,  la  diffé- 

S*<X> 

rentiation  de  l’égalité  m—l  y dx  par  rapport  au  temps,  sous  la 

d X 

condition  ts  = const.,  donnera 


ïmH~X  -dtdx=o' 


ou  bien,  en  substituant  finalement  à i sa  valeur  tirée  de  (i5)  et 
R.  1. 


IO 
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effectuant  l’intégration, 


07) 


\£dTtdx=a(: 


£1  l dj\ 

•ipoa  y dx  J 


Si,  pour  figurer  l’onde  par  une  courbe,  on  attribue  à chacun  de 
ses  éléments  dm,  d’abscisse  horizontale  actuelle  x et  de  longueur 
actuelle  dx  (prise  en  valeur  absolue),  une  ordonnée  verticale  re- 
présentative, égale  à la  condensation  correspondante  y,  et  que 
d’ailleurs  l’onde  ait  reçu  initialement  une  forme  simple  où  y,  à 
partir  de  zéro,  s’élève  peu  à peu,  depuis  la  tête  ou  le  front  de 
l’onde  jusqu’à  un  certain  sommet  ou  maximum  y*,  pour  décroître 
de  même  graduellement  de  y1  à zéro  entre  le  sommet  et  la  queue, 
les  deux  éléments  dm  de  même  hauteur  actuelle  y,  situés  l’un  en 
avant  et  l’autre  en  arrière  du  sommet,  auront  évidemment  des  vi- 
tesses de  propagation  inégales;  car,  dans  la  parenthèse  de  (17),  le 
terme  en  e*,  seul  différent  de  l’un  à l’autre,  sera  positif  pour  le  pre- 
mier, à pente  C~  négative,  et  négatif  pour  le  second,  à pente 

positive.  Il  est  clair  que,  par  l’effet  de  ces  célérités  inégales,  l’espa- 
cement des  éléments  dm  antérieurs  et  postérieurs  augmentera,  en- 
traînant, du  moins  en  moyenne,  une  augmentation  de  leurs  lon- 
gueurs dx  et  une  diminution  corrélative  de  leur  hauteur  y,  qui  est, 
dans  chacun,  inverse  de  la  longueur.  Ainsi,  grâce  à l’imperfection 
de  la  fluidité,  exprimée  dans  (17)  par  le  terme  en  e,  considéré, 
l’onde  s’allonge  et  s’aplatit,  comme  aurait  pu  le  faire  prévoir  la 
seconde  formule  (16),  où  l’intégrale  f y2  dx  = f y dm  ne  peut  dé- 
croître continuellement  sans  que  la  valeur  moyenne  de  y diminue. 

Mais  il  y a plus  de  complication  dans  la  manière  dont  se  dé- 
forment séparément  les  deux  parties  de  l’onde  qui  s’étendent, 
l’une,  de  la  tête  ou  du  front  au  sommet,  l’autre,  du  sommet  à la 

queue;  car,  suivant  le  plus  ou  moins  de  pente  qu’elles  offrent 

aux  divers  points  de  leur  longueur,  le  dernier  terme,  en  de  la 
parenthèse  (17),  peut  ou  non  l’emporter,  dans  des  mesures  va- 
riables aux  endroits  soit  bas,  soit  hauts,  sur  le  précédent,  simple- 
ment proportionnel  à y,  et  rendre  plus  prompte  à l’avant,  plus 
lente  à l’arrière,  tantôt  la  propagation  des  éléments  dm  de  faible- 
hauteur  y,  tantôt  celle  des  éléments  dm  plus  élevés. 


DES  ONDES  AÉRIENNES,  DANS  LES  TUYAUX  DE  CONDUITE.  n 
Lorsqu’on  admet  assez  de  petitesse  soit  dans  le  coefficient  e(, 
soit  dans  les  pentes  et  surtout,  par  suite,  dans  les  courbures 

pour  rendre  le  dernier  terme  de  (17)  ou  celui  de  (i5)  négligeables 
devant  les  précédents,  les  éléments  dm  qui  correspondent  aux 
plus  grandes  hauteurs  y sont  évidemment  les  plus  rapides;  et  le 
sommet  se  rapproche  sans  cesse  du  front  ou  s’éloigne  de  la  queue, 

en  adoucissant  de  plus  en  plus  les  pentes  ^ à l’arrière,  mais,  par 

contre,  en  les  accentuant  à l’avant  jusqu’à  ce  que  le  terme  en  e, 
y prenne  de  l’influence  et  les  empêche  de  devenir  abruptes,  ou  le 
sommet  d’être  à pic,  de  surplomber  presque. 

On  le  voit,  du  reste,  plus  simplement,  par  l’intégration  exacte 
de  l’équation  (i5)  où  le  terme  en  est  actuellement  négligèable, 
intégration  dès  lors  facile  qui  donne 

x — a|i  + m 1 yjt  ou  x — (a- h m 1 u}jt  = fond,  arbitr.  de  y ; 
ce  qui  oblige  à attribuer  une  vitesse  vers  les  x positifs, 


/ m + 1 \ m -+- 1 

a H y ) = a H u, 

\ a*/  2 ’ 

uniformément  croissante  avec  y (ou  avec  m),  à chaque  ordonnée  y, 
invariable  en  grandeur. 

Il  n’existe  donc  pas  de  forme  permanente  pour  l’onde,  quand 
on  néglige  les  frottements;  et  la  seule  qui  le  soit  lorsqu’on  en 
tient  compte  est  celle  qui  suppose  l’annulation  totale  de  l’énergie, 
ou  qui  a pour  équation  y = o,  forme  limite  n’exprimant  que  l’ex- 
tinction finale  du  phénomène.  Ainsi  une  intumescence  aérienne 
isolée  se  distingue  complètement,  sous  ce  rapport,  de  fonde 
liquide  solitaire  propagée  le  long  d’un  canal,  qui  comporte  et 
affecte  une  forme  permanente  bien  déterminée  pour  chacune  des 
valeurs  successives  de  son  énergie,  également  décroissante  sans 
cesse  (*  ). 


(J)  Il  n’en  serait  plus  tout  à fait  de  même  pour  une  intumescence  continue, 
c’est-à-dire  pour  une  onde  sans  queue,  où  la  hauteur  représentative  y,  nulle  à 
la  tête  x — 00  de  l’onde,  tendrait  peu  à peu  vers  une  cerlaine  valeur  2 c et  s’y 
maintiendrait  asymptotiquement  jusqu’à  x — — 00  (ce  qui  implique  une  énergie 
totale  infinie)',  car  il  suflirait,  pour  lui  conserver  sa  forme,  de  lui  donner  celle 
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VJI.  Mais  voyons  comment  le  frottement  de  la  paroi  du  tuyau 
et  sa  perméabilité  calorifique  modifieront  ces  lois,  en  annulant, 
sur  tout  le  contour  y.  des  sections  normales  <7,  la  vitesse  u du  fluide 
et  les  variations  relatives  t de  sa  température. 

Et,  d’abord,  la  petitesse  tant  du  coefficient  £ de  frottement  inté- 
rieur de  l’air,  que  de  sa  conductibilité  calorifique,  localiseront 
dans  une  couche  très  mince  les  perturbations  considérables  qu’en- 
traîne une  pareille  immobilisation  du  fluide  et  de  sa  température, 
de  manière  à laisser  subsister  très  sensiblement,  à l’intérieur  des 
sections,  la  forme  plane  des  ondes,  ainsi  que  la  longitudinal! té 
des  mouvements,  avec  leur  parité  du  centre  jusque  auprès  du  con- 
tour et  celles  de  la  condensation  y,  de  la  température  T,  enfin,  de 
la  pression  p.  Ce  ne  pourra  être,  du  moins,  que  dans  des  ondes 
très  longues  et  assez  loin  en  arrière  de  leur  tête,  c'est-à-dire  après 
une  action  prolongée  des  frottements  et  de  la  conductibilité  sur 
les  mêmes  particules  fluides,  que  les  perturbations  dont  il  s’agit 
auront  eu  le  temps  de  se  propager,  avec  une  intensité  notable, 
jusqu’à  des  distances  de  la  paroi  comparables  au  rayon  du  tuyau. 
Aussi  MM.  Violle  et  Vauthier  ont-ils  constaté  que  le  front  de  l’in- 
tumescence garde  indéfiniment  sa  forme  plane,  jusqu’aux  points 
les  plus  rapprochés  de  la  paroi  (om,o8)  où  ils  eussent  installé  un 
appareil  indicateur  de  l’arrivée  de  l’onde. 

On  le  conçoit,  du  reste,  d’autant  mieux,  que  la  forme  plane  s’y 
rétablirait  presque  instantanément  si  quelque  cause  accidentelle 
l’y  avait  altérée;  car  les  premières  particules,  centrales  ou  péri- 
phériques, atteintes  par  le  mouvement  sur  une  section  a-,  trans- 


qui  rend  l’expression  (17)  de  w indépendante  de  x ou  lu  même  d’un  bout  à 
l’autre  de  l’onde,  et  qui  assure,  par  conséquent,  avec  des  espacements  constants  dx 
entre  les  éléments  successifs  dz, j = y dx , la  permanence  de  leurs  hauteurs  y. 

La  parenthèse  de  (17)  s’y  réduisant  à 1 -+-  —~~~c  pour  x=  — oc  où  y=2 c^const., 

il  viendra,  comme  équation  du  profil  qui  définit  cette  forme, 


ou  bien 


/ N Po  CL  , 

— (m  + i)  - — dx, 


et,  en  intégrant, 


log 
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mettraient  aussitôt  ce  mouvement  au  devant  d’elles  près  des 
parois  non  moins  que  près  de  l’axe,  sur  toute  l’étendue  des  zones 
sphériques  de  rayon  croissant  les  ayant  pour  centres. 

Donc,  les  formules  (5)  à (11)  continueront  à subsister  dans 
l’intérieur  des  sections,  du  moins  à une  première  approximation. 

Cela  posé,  si,  considérant  un  élément  clf  du  contour  de  la  sec- 
tion quelconque  cr,  d’abscisse  x,  et,  en  même  temps,  une  portion 
contiguë  à aèy,  sensiblement  plane,  de  la  couche  superficielle  où 
sont  localisées  les  perturbations,  l’on  appelle  N,  dans  le  plan  de  a-, 
une  coordonnée  rectiligne  normale  à d’y,  introduite  pour  me- 
surer la  distance  à la  paroi  de  chaque  feuillet  de  cetle  portion  de 
couche,  la  vitesse  a et  la  condensation  y,  très  rapidement  variables 
à partir  de  la  paroi  ou  pour  les  très  petites  valeurs  positives  de  N, 
deviendront  sensiblement  indépendantes  de  cette  variable  dès 
qu’elle  atteindra  une  grandeur  appréciable;  et  elles  tendront  alors, 
à fort  peu  près,  vers  des  expressions  de  la  forme  (io),  savoir 
u = av(x  — at),  y=iv(x  — cit). 

Mais  formons,  au  même  degré  d’approximation,  les  équations 
indéfinies  du  mouvement  troublé,  pour  l’intérieur  d’un  pareil 
fragment  de  la  couche  superficielle. 


Les  composantes  transversales  e,  w de  la  vitesse,  milles, 
comme  u,  à la  paroi,  mais  nulles  de  plus  à l’intérieur,  sur  l’autre 
face  de  la  couche,  seront  négligeables  partout,  comparativement 
à w,  et,  de  même,  les  accélérations  latérales  e',  wf  disparaîtront  à 
côté  de  l’accélération  longitudinale  elle-même  réductible,  ici, 

à son  terme  principal  ou  linéaire  • Enfin  les  deux  dérivées  se- 
condes de  u soit  en  x et  y,  soit  en  x et  s,  pourront  être  négligées 
vis-à-vis  des  dérivées  secondes  de  u par  rapport  aux  coordonnées 

transversales  y ou  z,  vu  que  les  dérivées  ^ et  nulles  à l’in- 
térieur de  la  couche,  prennent  leurs  très  grandes  valeurs,  dans 
chaque  section  cr,  le  long  de  chemins  N extrêmement  petits,  mais, 
au  contraire,  varient  graduellement  suivant  le  sens  des  x.  Il  suit 
de  là  que,  dans  les  trois  premières  équations  indéfinies  (i)  du 
mouvement,  les  valeurs  (2)  des  pressions  N,  T,  différentiées  en  x , 


d2  u 


y , 3,  se  réduiront  aux  termes  en p et  aux  deux  termes  £ ^ ^ 2 dz2) 9 

, , f d2u  d2u\  , • 1 t , , . 

dont  la  somme  £ + ~JZ2)  Pourra  évidemment  s écrire  aussi  £ A2  u 
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et,  finalement,  Gomme  on  sera  d’ailleurs  libre  de  remplacer 

par  p0  les  dénominateurs  p des  premiers  membres,  il  viendra 

/ gx  ï_  dp  p d^u  _ du  i dp  _ i dp 

po  dx  p0  «?N2  dt  ’ p dy  ’ p dz  ° 

On  voit  donc  que,  même  près  du  contour  la  pression  p a ses 
dérivées  en  jk  et  z négligeables  à côté  de  sa  dérivée  en  x,  laquelle 
y est  de  l’ordre  de  sA2u.  La  pression  élastique  p , partout  propor- 
tionnelle à pT,  ou  égale  à po(i  -f-  y)(i  -b  t),  ou  enfin,  sensible- 
ment, à p0  (i  -f-  y -f-  t),  et  qui  d’ailleurs,  d’après  (5),  devient  à 
très  peu  près  p0  [i  H-  m <p(x — at )]  à l’intérieur,  sera  ainsi  con- 
stante sur  toute  l’étendue  des  sections,  même  dans  l'épaisseur  de 
la  couche  superficielle.  L’on  aura,  par  conséquent,  d’une  manière 
générale, 

(19)  Y + T = mcp(a?  — at). 


VIII.  Cela  posé,  pour  évaluer  d’abord  u à l’intérieur  de  la 
couche,  remplaçons  p , d’après  (6),  par  [1  _j_  cp(x  — at)\m  ou 
par  l’expression  sensiblement  égale  — h-o0a2<p(x — al ),  dans 
la  première  équation  (18);  et  cette  équation  deviendra 


, , , x £ d2  u du 

-a*,  (*-«*)  + ^3nî=** 

ou  bien,  en  observant  que  cp  ne  dépend  pas  de  N, 


(20) 


£ d2 


dt 


[u  — acpO  — a«)J=  “ — acp(x  — a*)]- 


Ainsi  écrite,  elle  a la  forme  de  l’équation  classique  des  tempé- 
ratures dans  un  mur,  et  contient  comme  fonction  inconnue 
l’excès  u — acp,  égal  à — acp  pour  N=  o,  mais  très  sensiblement 
évanouissant  dès  que  N acquiert  des  valeurs  positives  perceptibles. 
Les  valeurs  notables  de  la  fonction  u — ■ acp,  réalisées  toutes  dans  le 
voisinage  de  la  paroi  ou  de  N = o,  seront  donc  à fort  peu  près  celles 
de  la  température  près  de  la  face  chauffée  d’un  mur  très  peu  con- 
ducteur, dont  la  masse  intérieure  aurait  anciennement  été  portée  et 
resterait,  pendant  un  temps  presque  indéfini,  à la  température  zéro, 
tandis  que  sa  couche  superficielle  subirait,  sur  sa  face  externe,  -la 
suite  des  températures  données  — av(x  — at),  sensiblement  pa- 
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reilles  à un  même  instant  dans  tout  le  voisinage  de  l’élément  de 
contour  dy.  L’intégrale  (connue)  de  (20)  qui  donne  u — acp  dans 
ces  conditions,  est 


(21) 


u — a cp(a?  — at)  — 


^x  — at  - 1- 


«PoN«\ 
2 sa2  / 


_ a* 

e 2 da, 


comme  on  peut  le  voir  par  lesnos346*  et  466*  de  mon  Cours  d’A- 
rtalyse  infinitésimale  (t.  II,  p.  178*  et  468*),  où  il  est  démontré 
aussi  que  la  dérivée,  par  rapport  à N,  de  u ou  du  second  membre 
de  ( 2 1),  égale 


a2  p 0 

E 


1 ru 


x — at  -+- 


ap0a2 
2 £ 


_ Ai 
e 2a*  da.. 


Bornons-nous  à en  déduire  le  frottement  § qu’exerce,  dans  le 
sens  des  x négatifs,  l’unité  d’aire  de  la  paroi,  sur  la  couche  gazeuse 

contiguë,  frottement  dont  l’expression  est  Il  faudra  donc 


prendre  cette  expression  à la  limite  N — o,  et  il  viendra,  en  chan- 
geant finalement  la  variable  d’intégration  de  manière  à poser 


(22) 


x — at 


ap  o«2 
2 £ 


i'#  = — 

= — 2ap0i/—  f v'(x — at->r^)d^. 
V 'TCpo  J0 


da. 


IX.  Passons  maintenant  au  calcul  analogue  des  variations  rela- 
tives t de  la  température,  sur  les  divers  feuillets  de  la  même  portion 
de  couche  gazeuse  contiguë  à dy , en  vue  d’évaluer  les  conden- 
sations y,  reliées  à t par  l’équation  (19). 

A la  paroi,  c’est-à-dire  pour  N = o,  l’on  aura 

t = o,  ou  y = m y (x — at) 

et,  à l’intérieur,  c’est-à-dire  dès  que  N acquerra  des  valeurs  po- 
sitives perceptibles,  y tendra  sensiblement  vers  <p(#  — at ),  ou  t, 
d’après  (19),  vers  (m  — i)cp(.r — at)  • L’excès  y — cp(#  — at) 
variera  donc  très  vite,  comme  u , depuis  une  valeur  notable  censée 
donnée,  qui  est  ici  (m  — 1)  cp(.r  — at),  jusqu’à  zéro,  quand  la  va- 
riable N,  d’abord  nulle,  passera  elle-même  par  toutes  les  très 
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petites  valeurs  positives.  Or  la  relation  (4)  [n°  III]  fournit  aisé- 
ment l’équation  indéfinie  nécessaire  pour  achever  de  le  déter- 
miner dans  cet  intervalle. 

Appliqués  à un  feuillet,  d’épaisseur  c/N  et  de  masse  M=p  cTN  par 
unité  d’aire,  du  fragment  considéré  (contigu  à dy)  de  couche  su- 
perficielle, les  premier  et  quatrième  membres  de  celte  relation  (4), 
divisés  par  dt , donneront,  en  y remplaçant  dT  par  T0  d'z,  c/p  par 
p0  et  puis,  à fort  peu  près,  p,  T par  p0,  T0, 

(23>  §=CT.Po«iN  [*-(«-: ,)$]• 


D’ailleurs,  la  quantité 


dQ 

dt 


de  chaleur,  reçue  par  le  feuillet  durant 


l’unité  de  temps,  ne  sera  plus  négligeable  comme  pour  les  parti- 
cules intérieures,  malgré  l’extrême  petitesse  du  coefficient  k de 
conductibilité  de  l’air;  car,  dans  l’expression  des  deux  flux  calo- 
rifiques qui  traversent  les  deux  bases  du  feuillet  (suivant  le  sens 
inverse  de  celui  des  N croissants),  cette  petitesse  du  facteur  k 

sera  compensée  par  l’énorme  grandeur  de  la  dérivée  — ou 

-b  “^■7  e/N  le  multipliant.  L’on  aura,  par  conséquent,  toujours 

pour  l’unité  d’aire  du  feuillet,  en  prenant  la  différence  de  ces  deux 
flux  après  avoir  donné  à k , sauf  erreur  négligeable,  sa  valeur  k0 
aux  température  et  densité  primitives  T0,  p0, 


dQ 

dt 


L’équation  (23),  divisée,  CT0 p0  û?N,  deviendra 


(24) 


ko  d*  t _ çfc  _ v rfy 

p0CdN*“<fc  1 dt  * 


Remplaçons-y  enfin  t par  sa  valeur  tirée  de  (19)  et  nous  aurons, 
pour  régir  la  différence  y — — at)  dans  toute  l’épaisseur  de 

la  couche  superficielle,  l’équation  indéfinie  cherchée,  analogue 

à (20), 

(25)  *)]  = ^55Ê  [*-*(*-“*>]• 


Intégrons  celle-ci  de  la  même  manière  que  plus  haut  l’équation 
(20),  en  observant  que  la  valeur  de  la  fonction  inconnue  pour 
N = o est  maintenant  (m  — 1)  <p(.r  — at). 
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La  formule  de  cette  fonction  sera,  pareillement  à (21), 


(■26)  y — y(x  — at)  = (m  — 1) 


1 Ajf 


al 


«rap0CN2 ^ 
2 An  a2 


a * 
2 


dN 


Nous  aurons,  ci-après,  à en  déduire  la  valeur  de  l’intégrale 

d 
'0 

qui,  d'après  (a5),  est 

,,,  x-,  /*\ 


/ a tT -?(*-“)] 


A-o 

po  C J0  dN* 


dN  = - ) , 

mp  0C  \aN/o 


A, 


et  se  réduit,  par  suite,  au  produit  de  — — par  la  dérivée  du 

second  membre  de  (26)  prise  à la  limite  N — o.  En  évaluant  cette 
dérivée  comme,  plus  haut,  celle  du  second  membre  de  (21),  et 

puis  posant  finalement  a = S 4 / — x » il  viendra 
1 J 1 y amp0  C 


(27) 


. t A /**  ,/  amp0Ca2\ 

= cpp-a<+  J «fa 


X.  Maintenant  que  nous  avons  exprimé,  en  fonction  des  valeurs 
de  Met  de  y à l’ intérieur  de  la  masse  gazeuse,  le  frottement  exercé 
d sm  surface  et  les  conden sations  spéciales  qu’elle  y éprouve,  nous 
pourrons  former,  à une  seconde  et  suffisante  approximation,  les 
équations  de  mouvement  de  tonte  une  tranche  fluide,  comprise, 
à l’époque  t , entre  les  deux  sections  normales  <7  du  tuyau  qui  ont 
les  abscisses  respectives  x et  x -f-  dx. 

Appliquons  d’abord  à cette  tranche  le  principe  des  quantités  de 
mouvement,  suivant  le  sens  des  x , en  écrivant  que  la  somme 

dx  / M'(p  dn ) des  produits  des  masses  élémentaires  p d?  dx  qui  la 
d a 

composent,  par  leurs  accélérations  longitudinales  u\  égale  celle 
des  actions  extérieures  exercées  du  dehors  sur  la  tranche,  suivant  le 
même  axe  des  x \ ce  qui  comprend,  d’une  part,  la  somme  algébrique 
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d r 

dx  ~r  I N xd<r  des  actions  normales  sollicitant  ses  deux  bases 
dx  J « 


et,  d’autre  part,  le  frottement  total  changé  de  signe,  — dx  I -f  dy, 
que  subit  sa  tranche. 

Dans  les  premières  de  ces  forces,  la  partie  principale  de  N*  sera 
la  pression  élastique  changée  de  signe,  — p,  constante  sur  toute 
l’étendue  cr  de  chaque  section,  et  exprimable  par  la  formule  (6) 
[n°  III]  où  nous  pourrons  supposer  la  valeur  de  y relative,  par 
exemple,  au  centre  des  sections.  Quant  à l’autre  partie  de  N*., 
due  à l’imperfection  de  la  fluidité,  et  tout  entière  de  seconde 
approximation,  elle  pourra  être  réduite  très  sensiblement  à sa 


valeur  dans  l’intérieur  de  la  section,  savoir  à £ 


du 


> '^0,ia 


d.a<o(x — at ) , 

61 — ‘s — = E*a<?> 

quantité  dont  elle  ne  diffère  dans  un  rapport  notable  que  sur  une 
étroite  zone  périphérique  de  <7,  relativement  négligeable.  Enfin, 
comme  la  valeur  (22)  de  § sera  la  même  pour  tous  les  éléments  dy 
du  contour  d’une  section,  le  second  membre  de  l’équation  s’écrira 

( =-a<*r|^p<la2(i  + Y)'n-1^  — 

On  le  simplifie  encore,  en  observant  que  p0«2(  1 H- y)m~  1 n’excède 

p°a2^que  d’un  terme  de  deuxième  approximation,  réductible 

à p oa^Çm  — i)y  ou  même  à p 0a2(m  — 

Quant  au  premier  membre  dx  / u'pdv,  ou  pQdx  / (ur-hyu')d<r, 

d(j  J a 

la  partie  linéaire  et  principale  de  u'-\-yu'  y est  celle  de  u\  savoir 

et  son  autre  partie,  de  deuxième  approximation,  peut,  comme  plus 
haut  celle  de  N*,  n’être  prise  que  dans  l’intérieur  de  <7,  où  elle  se 

réduit  sensiblement  à u +Y~^r  et  même  à a2cpcp'-{-<px( — a1  cp')=o. 

Appelons,  d’ailleurs,  U la  vitesse  moyenne  exacte  à travers  la 
section  <7  (moyenne  très  peu  différente,  évidemment,  des  valeurs 
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de  u à l’intérieur,  ou  même  de  acp),  de  manière  à pouvoir  remplacer 

r du  t d r j d. U a d\ J t ,,  , . , , , 

J dv  ou  J u a<7  par  = a-  et  1 équation  cherchée, 

divisée  par  p0  vdx,  sera  enfin,  grâce  à la  dernière  expression  (22) 
de 


(28) 


dû 

~dt 


9 d'[  d r 

° dx  dx  L 


m — 1 


a2  y(x  — at)'2  — — a o'(x  — at ) 

Po  1 J 

f ~ at  + P2)  = O* 


La  seconde  relation  nécessaire  pour  déterminer  les  variations 
successives  des  deux  fonctions  inconnues  U et  y est  fournie  par 
l’équation  dite  de  continuité , qu’il  suffira  d’appliquer  à l’espace 
compris,  dans  le  tuyau,  entre  les  deux  sections  fixes  d’abscisses  x 


et  x + dx.  La  masse  gazeuse  dx  / p d<7  ou  p0dx  I (î-f-y ) â?<7, 

d (y  d (y 

remplissant  cet  espace  à l’époque  t , s’accroît  dans  l’unité  de  temps 
de  l’excédent,  — dx—J'pudcr,  de  la  matière  J*  p u d<j  entrée 


par  la  première  section  sur  celle,  d’expression  analogue,  qui  est 
en  même  temps  sortie  par  laseconde.  Et  comme  p w = p„(w  + wy), 
il  vient  ainsi,  après  division  par  p0  dx , 


(29>  + 

Or,  dans  le  calcul  du  troisième  membre,  le  dernier  terme,  de 
seconde  approximation,  peut  être  évalué  en  négligeant  les  per- 
turbations localisées  sur  le  contour  et  en  remplaçant  même  y,  u 
par  <p(x  — at ),  az>(x  — at).  D’autre  part,  dans  le  calcul  du  pre- 
mier membre,  l’aire  totale  <7  comprendra  : i°,  l’élroite  zone  péri- 
phérique, où  la  condensation  y pourra  être  remplacée  par  sa  va- 
leur de  première  approximation,  égale  à celle  de  y au  centre, 
augmentée  de  l’excédent  de  y sur  (p(x  — at)  qui  exprime  sensible- 
ment cette  valeur  dey  au  centre,  et,  20,  la  partie  intérieure,  prin- 
cipale, où  y,  devenue  indépendante  de  la  distance  N au  contour, 
est  à très  peu  près  la  valeur  de  y au  centre,  que  nous  désignons  sim- 
plement par  y.  Le  premier  membre  de  (29)  se  dédouble  ainsi  en 
deux  termes,  dont  le  plus  simple,  qui  à lui  seul  l’exprimerait  sans 

les  perturbations  produites  au  contour,  est^-  C d(T  = iJ^~y  tandis 
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que  le  second,  composé  d’éléments  de  la  forme 

^ Eï  — ?(*-«*)] 


sera  évaluable  en  divisant  la  zone  périphérique  en  rectangles  élé- 
mentaires dN  par  des  normales  et  des  parallèles  au  contour. 

Gomme,  d’ailleurs,  pour  la  bande  qui  correspondra  à chaque  élé- 
ment d du  contour,  l’intégration  par  rapport  à N pourra  se  faire 
de  zéro  à l’infini  [à  cause  de  l’annulation  asymptotique  de 
y — — at)  dès  que  N devient  sensible],  la  formule  (2^ ) four- 

nira la  valeur  de  ce  second  terme,  et  les  deux  membres  extrêmes 
de  (29),  divisés  par  a-,  donneront  enfin  l’équation  cherchée,  qu’il 
faudra  joindre  à (a 8), 


(3o) 


I dv  d\ J d. o(cr  — aty 

j dï+dï+a  ‘~~cho 

I — 2 ~(m — 01/  — / v’(x  — a/ -t- 82)  = o. 

\ <x  y T:mp0  LJ  ‘ 


Tandis  que  les  deux  équations  (28)  et  (3o)  du  problème 
contiennent  la  vitesse  moyenne  U sur  une  section  a-  et  non  la 
vitesse  à l’intérieur,  c’est  la  condensation  y produite  à l’intérieur 
que  la  formule  (6)  de  la  pression  p nous  a fait  introduire  dans  la 
première,  (28),  et  que  nous  avons  employée  aussi,  par  suite, 
dans  (3o),  à la  place  de  la  valeur  moyenne  de  y sur  toute  la  sec- 
tion a-  d’abscisse  x , qui  aurait  figuré  plus  naturellement  dans  le 
premier  membre  de  (3o),  en  le  débarrassant  de  son  dernier  terme, 
le  plus  compliqué.  Si  nous  appelons  T cette  moyenne  des  valeurs 
de  y sur  toute  une  section,  le  produit,  par  <7,  de  son  excédent  sur 
la  valeur  de  y à l’intérieur,  égalera  évidemment  l’intégrale  de 
l’expression  [y  — v(x  — at )]  dN  dj^  prise  pour  tous  les  éléments 
d?  — dN  dj^  de  la  zone  périphérique,  intégrale  qui,  obtenue  en 
se  servant  de  la  formule  (27)  intégrée  de  t = — 00  à t=t,  est 
évidemment 


Ou  aura  donc,  pour  remplacer,  au  besoin,  y par  Y dans  les  équa- 
tions du  problème,  la  relation 
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XI.  On  remarquera  que,  — at ) eta<f(x — at)  exprimant 

y et  U dans  tous  les  termes  de  deuxième  approximation,  les  pre- 
miers membres  des  deux  équations  (28),  (3o)  ne  diffèrent  de  ceux 
des  précédentes  (12)  [n°  IV]  que  par  l’addition,  dans  chacune, 
du  dernier  terme,  où  figure  une  même  intégrale  définie.  Ces  deux 
termes  sont  donc  les  seuls  par  lesquels  se  traduise  l’influence  tant 
retardatrice  qu’extinctrice  de  la  paroi.  Et  ils  se  réduisent  même 
à un  lorsque,  après  avoir  respectivement  difîerentié  en  x et  t , 
puis  multiplié  par  — 1 et  + 1 , ces  deux  équations  (28),  (3o),  on 
les  ajoute  pour  éliminer  U.  11  vient  alors  en  y,  comme  dans  le  cas 
idéal  d’une  paroi  sans  frottement  et  imperméable  à la  chaleur, 
l’équation  (i3)  [n°  IV],  mais  où  K,  quantité  encore  nulle  à la 
tête  de  l’onde  (pour  x = 00),  a l’expression 


K 


(32) 


d_ 

dx 


d v(x  — at) 


[m  H-  t a . 

‘<22cp(a?  — at )2 — ~ , 

2 1 p o dx 

— 2 i*  £a  J*  v(x  — at  -h  P2)  d$ j , 


avec  un  nouveau  coefficient  p défini  par  la  relation 


(33) 


-4-  ( m — 1 ) 


Dans  le  cas  de  l’air  à la  température  ordinaire  et  sous  la  pres- 
sion atmosphérique,  on  a,  d’après  MM.  Violle  et  Vauthier  (p.  20 
de  leur  Mémoire),  en  prenant  le  mètre  et  la  seconde  pour  unités 
de  longueur  et  de  temps, 


0,0039, 


0,0057  ; 


et  l’on  sait  d’ailleurs  que  m = 1 , 4 ou  que  \m  — 1 » i832.  Il  en  ré- 
sulte 

p.  = o,oo58, 


valeur  assez  forte,  si  le  tuyau  n’est  pas  d’un  diamètre  de  l’ordre 
de  a ou  que  le  rapport  ^ s°it  considérable,  pour  rendre  le  der- 
nier terme  de  la  parenthèse  de  (32)  très  supérieur  au  précédent, 
où  se  trouvera  comparable  à ~ et  atteindra  seulement  quelques 
unités  du  cinquième  ordre  décimal. 
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L’équation  en  y du  premier  ordre  par  rapport  au  temps  sera 
donc  encore  (i4),  c’est-à-dire,  au  lieu  de  (i5), 


dy  dy  û?  |m  + i . £1  dv(x 

dt  + adï  + d*  I — «,(»-«*)»-  — -J- 


at) 


(34) 


dx 


<p(# — at  -F  P2)  = o. 


On  pourra  y remplacer,  au  besoin,  y par  T ; car  la  Formule  (3i) 
montre  que  la  différence  T — y est  une  fonction  de  x — at  ou  a sa 
dérivée  en  t égale  au  produit,  par  — a,  de  sa  dérivée  en  x. 

D’ailleurs,  la  remarque  précédente  permettra  généralement  de 
supprimer,  dans  la  quantité  entre  crochets  de  (34),  le  terme  en  . 
On  pourra  même  y effacer  le  précédent,  affecté  du  carré  cp2  de  la 
condensation  y,  sauf  au  commencement  de  la  propagation  où 
l’onde  n’est  pas  encore  affaiblie  et  où,  comme  nous  avons  vu  au 
n°  VI,  l’effet  de  ce  terme  est  de  tendre  à rapprocher  le  sommet  de 
l’intumescence  de  son  extrémité  avant. 


XII.  Ici,  la  masse  de  l’onde,  représentée  proportionnellement 
par  f y dx  ou  mieux  par  fV  dx , n’a  plus  la  même  importance  que 
dans  le  cas  d’intumescences  rigoureusement  planes,  parce  que 
l’action  des  parois  ne  tarde  pas  à rendre  inefficace  une  partie  de 
plus  en  plus  notable  de  cette  masse  totale,  en  la  disséminant,  à 
l'arrière  du  gros  de  l’onde,  sur  une  grande  étendue  où  la  con- 
densation y est  insensible  (*)  et  où  ne  se  conserve  plus,  dans  le 
fluide  presque  en  repos,  la  parité  des  mouvements  supposée  par 
nos  formules  à l’intérieur  de  chaque  section  normale. 

Mais  à cette  portion  disséminée  de  la  masse  de  l’onde,  il  ne 
correspond  qu’une  énergie  totale  insignifiante  : car,  vu  la  double 
proportionnalité  (sous  volume  égal)  des  éléments  de  masse  à y et 
des  éléments  d’énergie  à a1  y2,  ceux-ci  ne  sont  en  tout,  pour  la 
partie  en  question,  que  de  l’ordre  de  sa  masse  finie  multipliée  par 
la  valeur  moyenne,  négligeable , qu’y  présente  y. 


( 1 ) L’inévitable  longueur  de  cette  queue  se  reconnaît,  sur  l’équation  ( 34  ) multi- 
pliée par  dx  et  intégrée  depuis  l’abscisse  actuelle^  d’un  de  ses  points  jusqu’à  x = ac, 

à ce  que  la  dérivée  en  t de  /y  dx  qui  en  résulte,  — p.  — i - / y(x. — a£-b32)<73. 

y * j o • 

est  négative  (pour  tp  > o),  et  non  pas  nulle,  tant  que  x0  ne  devient  pas  — oc. 
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Donc  presque  toute  l’énergie  tant  actuelle  que  potentielle  du 
mouvement,  encore  subsistante,  réside  dans  le  gros  de  l’onde  non 
disséminé  auquel  s’appliquent  nos  formules,  et  a son  expression 
sensiblement  proportionnelle  à l’intégrale  /y2  dx,  calculée  entre 
les  limites  # = ±00  pour  la  fonction  y que  régit  l’équation  (34). 
Or  celle-ci,  multipliée  par  2y dx,  puis  intégrée  terme  à terme, 
entre  les  limites  dz  oo,  comme  on  a fait  pour  obtenir  la  seconde 
formule  ( 1 6 ),  donne,  en  réduisant  d’ailleurs  y à <p(#  — at)  dans 
les  termes  de  deuxième  approximation, 

d r"  s.  z*00 

— I y2dx  = / cp'(a? — at)2  dx 

dtJ_J  Po 

^ f ^^x~at^^x~at+^dx’ 

Telle  est  la  relation  qui  remplacera  la  seconde  (16)  pour  éva- 
luer le  décroissement  incessant  d’énergie  de  l’onde.  Malheureu- 
sement le  second  membre,  même  simplifié  par  la  suppression  du 
terme  en  permise  d’ordinaire  (comme  on  a vu),  ne  paraît  pas 
moins  difficile  à interpréter  d’une  manière  générale  que  le  dernier 
terme  entre  parenthèses  de  (34).  Aussi  nous  bornerons-nous  à voir 
ce  que  donnent  l’une  ou  l’autre  de  ces  équations,  mais  surtout  la 
première,  plus  compréhensive,  (34),  dans  des  cas  où  s’effectuent 
simplement  les  intégrations  qu’elles  impliquent. 

XII 1.  Le  premier  de  ces  cas  est  celui  de  mouvements  devenus 
depuis  longtemps  pendulaires,  avec  un  certain  nombre  constant,  /r, 
d’oscillations  complètes  par  seconde.  On  conçoit  qu’alors  la  forme 
plane  des  ondes  se  conserve  indéfiniment,  à l’intérieur,  malgré 
l’influence  des  parois;  car  celle-ci,  s’exerçant  dans  des  sens  alter- 
nativement opposés,  comme  ceux  des  vitesses  et  des  variations 
de  densité  qui  la  provoquent,  ne  peut  s’accroître  d’une  période 
d’oscillation  à l’autre.  L’équation  (34)  reste  donc  applicable  même 
quand  la  tête  ou  le  front  de  l’onde  s’est  propagé,  pour  ainsi  dire, 
jusqu’à  une  distance  infinie  x de  l’origine  des  coordonnées,  et 
quand  l’on  a,  par  suite,  sensiblement,  à un  facteur  constant  près, 
pour  toutes  les  abscisses  x de  grandeur  modérée,  en  adoptant 
comme  origine  des  temps  une  époque  où  s’observe  aux  points  x = o 
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un  maximum  de  condensation, 

2mz(x  — at  ) 


(36) 


cp(a?  — at)  = cos 


Alors  la  parenthèse  de  (34),  réduite  d’ailleurs  à son  terme  en  p, 
devient 


nz(x — at)  2/i7r(32 


*Wli  [™- “• 


. 2mz(x — at) 


- -!£=*[< 
i \J  n TT  a L 


. 2/l7Tp2‘] 

un  — t/p 

«J 


2iiiz(x  — at)  . txmz{x  — at)~ 


L’équation  (34)  a donc,  pour  son  troisième  terme, 

y \ . 2itk(x — at)  ‘inizix  — at)~\ 

— — sin -h  cos 

•2  a L a a 


ou  encore 

fj.  y/ Ai  TC 


ti  71  x r t 

2 a L 2 /l  7T 

et  enfin,  plus  simplement, 


2mz(x — at)  2 /nr(r  — «O 

■y-  cos h COS 

tc  cit  a a 


iW^  Z A,  1 ch 

2 a \ 1 2 /i7c  dt 


Cette  équation  (34)  peut,  dès  lors,  s’écrire 


(37) 


X.\  d1 


\j  niz  a 


[jl  \J  n 7T  y 


a -7-  h-  •— - 
ctx  2 a 


et  son  intégrale  générale  est  de  la  forme 


(38) 


ui/niT  7 

avec  m = -,  o>  = 

2 a a 


f = e-'»*  f(x  — a)t), 
a 


i -b 


Z 

4 \//iTr  a 


Z 

4/nx  » 


27C  R 


Si  l’on  remplace  ^ par  ou  par  ~ (en  supposant  ainsi  les 

tuyaux,  circulaires  et  de  rayon  R),  dans  ces  deux  formules  de  m et 
de  (o,  elles  deviennent  précisément  celles  que  KirchhofF,  et  anté- 
rieurement M.  Helmholtz  ^sauf  la  réduction  de  p à y/ avaient 
obtenues  par  de  tout  autres  méthodes.  Le  coefficient  d}extinc- 
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tion  m et  le  ralentissement  ou  retard  a — ou  sont  donc  bien,  le 
premier,  en  raison  directe  et,  le  second,  en  raison  inverse  de  la 
racine  carrée  du  nombre  n d’oscillations  par  seconde,  c’est  à-dire 
en  raisons  respectivement  inverse  et  directe  des  racines  carrées  de 
la  période  de  vibration  ou  de  la  longueur  d’onde.  Le  coefficient 
d’extinction,  en  particulier,  croît  quand  la  période  diminue;  et 
ces  ondes  aériennes  s’éteignent,  comme  les  ondes  liquides,  d’autant 
plus  vite  qu’elles  sont  plus  courtes  (1). 

XIV.  Mais  notre  but,  ici,  est  de  considérer  spécialement  les 
intumescences  limitées,  pour  voir,  par  exemple,  comment  s’abaisse 
peu  à peu  leur  sommet  (ou  comment  y décroît  la  condensation  y 
maxima)  et  aussi  de  combien  la  vitesse  de  propagation  de  ce 
sommet,  alors  que  les  frottements  l’ont  déjà  fort  abaissé,  se  trouve 
excédée  par  la  vitesse  a du  son  à l’air  libre. 

Intégrons,  à cet  effet,  l’équation  (34),  à partir  d’une  époque 
quelconque  tQ,  pour  un  temps  t — t0  assez  modéré,  de  manière 
que  la  solution  de  première  approximation  y = <p(#  — at)  y reste 
effectivement  approchée  avec  la  forme  de  la  fonction  <p  à l’é- 
poque c’est-à-dire  en  admettant  que  l’expression  exacte  de  y à 
cette  époque  soit  <p(# — at0).  Alors  la  relation  (34),  d’ailleurs 
réduite  à l’équation  du  premier  ordre 

(39)  M+aÈ  = ^\/îf0 

aura  pour  son  intégrale  générale 

(40)  Y = cp(|)-t-(^ — + où  il  est  posé  \ —x  — at. 

La  dérivée  en  x de  y,  nulle  au  sommet  de  l’onde,  sera  donc 

(40  è + 

Appelons  la  valeur  de  i*  qui,  à l’époque  t0 , correspond  au 
sommet  considéré  et,  par  suite,  au  maximum  de  cp(£). 


(*)  II  ne  serait  pas  difficile  de  reconnaître  la  parfaite  concordance  du  coeffi- 
cient d'extinction  m donné  par  la  seconde  formule  (38),  avec  celui  d’une  houle 
à vagues  longues  (comparativement  à la  profondeur  d’eau)  qui  se  déduit  de  la 
formule  (28)  du  Mémoire  cité  du  Journal  de  Mathématiques  (t.  IV,  p.  355; 
1878). 
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Ce  sommet  se  déplaçant  avec  une  vitesse  presque  égale  à <z, 
l’excédent  £ de  son  abscisse  x sur  at  ne  dépassera  , au  bout 
du  temps  t — t0 , que  d’une  quantité  du  premier  ordre  de  peti- 
tesse, et,  vu  l’annulation  de  cp' ( Ç , on  aura  très  sensiblement 
<p7(5)  = cp^Si  )(x  — at  — ).  D’autre  part,  le  terme  <p(£)  de  (4o) 

ne  différera  du  maximum  voisin  <p(Ç,)  que  d’une  quantité  négli- 
geable du  second  ordre  et  enfin,  dans  les  petites  parties  en  pi, 
de  (4o)  et  (40?  5 se  trouvera,  de  même,  réductible  à Appli- 
quées au  sommet  de  l’onde,  les  deux  équations  (4°)  et  ( 4 1 ) de- 
viennent donc 


(42)  < 0 ^ 

| o-?"(£i)(^  — at—  îi)  + (*'—  ?"(?i  + P2)^?- 

Différentions  la  première  par  rapport  à t , après  l’avoir  divisée 
parla  valeur  même,  o(Ç4  ),  de  y à l’origine  t0  de  l’époque  considérée. 
Il  viendra,  pour  régir  le  décroissement,  à cette  époque,  du  maxi- 
mum, l’équation  différentielle 


(43) 


1 

T 


dp. 


Son  intégrale,  si  l’on  appelle  am  le  second  membre,  constant  au 
moins  pendant  un  temps  modéré,  est  = const.,  relation 

exprimant  que  la  condensation  maxima  décroît  comme  l’exponen- 
tielle e~ma<'t~t 0),  pendant  qu’elle  se  propage  sur  un  parcours  sen- 
siblement égal  à af  — <£0).  Donc  m n’est  pas  autre  chose  que  le 
coefficient  cT extinction  cherché*,  et  si,  d’autre  part,  l’on  différentie 
également  en  t la  seconde  équation  (42)  divisée  par  <Z3"(I*,),  en 
appelant  w la  dérivée  de  x ou  vitesse  de  propagation  du  sommet, 
les  deux  formules  (pour  le  sommet  de  l’onde)  du  coefficient  actuel 

d' extinction  et  du  retard  relatif  actuel  a M » seront 


(44) 


y/  Cl  TC 


x r 

<7  L 


— ?'(Êi 


32) 


?(Êi) 


dp, 


a~  = X C*  ?"(£i-+-  P2)  d% 
\ a \J  cm  5 Jq  ? (ît  ) 


XV.  Dans  le  cas,  par  exemple,  d’ondes  périodiques,  où  la  fonc- 
tion c p a sensiblement  l’expression  (36),  on  peut  prendre  = o et 
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les  deux  formules  (44)  donnent  presque  immédiatement  les  va- 
leurs (38)  de  m et  de  to. 

Mais  attribuons  plutôt  à la  fonction  cp(£)  une  forme  qui,  tout  en 
rendant  facile  le  calcul  des  seconds  membres  de  (44 ) ? puisse  con- 
venir à une  intumescence  isolée,  ou  s’annule  asymptotiquement 
aux  deux  limites  ^ = dzoo  sans  présenter,  dans  l’intervalle,  aucune 
discontinuité.  La  plus  simple,  en  y faisant  abstraction  du  facteur 
constant  qui  exprime  la  condensation  maxima  initiale  (relative  à 
l’époque  t0 ),  me  paraît  être  celle-ci,  symétrique  de  part  et  d’autre 
de  la  valeur  de  sa  variable,  = o,  pour  laquelle  se  produit  son 
unique  maximum  i, 

(45)  = 

c y désigne  une  ligne  donnée,  en  rapport  avec  la  longueur  appa- 
rente de  l’onde  à l’époque  considérée  t0.  Si,  introduisant  une  va- 
riable auxiliaire  r\  qui,  de  la  tête  £=oo  à la  queue  ç=  — oo  de 

1 ? * A 1 / ' ^ 

I intumescence,  croîtra  de  zéro  a -»  nous  posons 

(46)  £ = ccot27),  on  aura  <p(jj)  = sin2  2r(,  ®(î)^E= — 2C  dr,; 

et  l’intégrale  f — atù)‘1dx  = /<p(£)2  d\,  représentative  de  l’é- 
nergie initiale,  sera  simplement  [vu  la  valeur  moyenne  ^ de  cp(£)  ou 
de  sin2 27]]  la  moitié  de  l’intégrale  / cp(ç)<f£  ===  — 2cf  dr{  = tcc,  re- 
présentative de  la  masse  sensible  ou  concentrée  de  l’intumescence. 

II  viendra  donc 

(46  bis)  f cp(£)2  de,  = - ne. 

Or,  pour  la  moitié  convexe  d’une  onde  sinusoïdale  où  <p  aurait 
l’expression  (36),  la  valeur  moyenne  de  cp2  serait  justement  et 
l’intégrale  f cp(ç)2  d\  y vaudrait  d\,  c’est-à-dire  la  moitié  de  la 
longueur  de  la  demi-onde  considérée.  Appelons  L cette  longueur, 
quotient  de  la  vitesse  a du  son  par  le  nombre  m des  oscillations 
simples  durant  une  seconde,  et  l’intumescence  exprimée  par  (43) 
sera  l’équivalent,  en  énergie,  de  la  partie  convexe,  de  même  hau- 
teur, d’une  onde  sinusoïdale,  qui  aurait  sa  longueur  totale  2L,  ou 
l’inverse  n de  sa  durée  totale,  liés  au  paramètre  c par  les  formules 

L _ a 

7T 


(47) 


2 /ITT 


(48) 
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Cela  posé,  évaluons  de  zéro  à go  l’intégrale  f cp  ( Ç — f-  en 

calculant  d’abord  cp ( Ç —1—  ) par  le  second  membre  de  (45),  et 

puis  l’intégrale  indéfinie  par  la  méthode  classique  (après  avoir 
substitué  c cot  2'r\  à Ç).  Il  vient,  comme  une  différentiation  immé- 
diate permet  d’ailleurs  de  le  vérifier, 


r P 

/ ?(É 

•■'O 

= \/ 


P2)<*P 


/(Mi 


, . , . /vrv-  . sJ C tang7))2+  c cot  7) 

- CO  S3  Y)  SUIT)  log4  / )■■■  ° — - 

2 V (Pv2  — \' c tangTj )2-b  c cotrj 


V\ 


sin37)  cos Tj  I arc  tang 


+ y/c  tan 

\/ C COtT) 


g 7)  VJ 

— -+-  arc  tang 


S v/2  — \/c  tangrj 


sJ c cot  7) 


expression  finie  et  bien  continue  pour  toutes  les  valeurs  de  [3, 
quel  que  soit  l’arc  7 ) entre  zéro  et  11  suffit  donc  d’v  faire  ^ = go, 
et  l’on  a 


(49) 


j ' cp(£  + (32)  6^(3  = sin 


7)  cos  - 7). 


Différentions  deux  fois  par  rapport  à ç,  en  nous  souvenant  que 

la  dérivée  de  r\  en  H est  — 'îlï)  0u  — -sin2r,  cos2tï;  et,  si  nous 

2 c c 

faisons  finalement  £ = = o,  7j  = — > siny^  = cosr,  = il  viendra 

4 y 2 

(5o)  r<f'at+p).dp= — p) ap =--??=■ 

4 \‘2C  do  8cy2C 

Comme,  d’ailleurs,  on  trouve  aisément  que  cp(£)  vérifie,  en 
vertu  de  (45),  l’équation  différentielle  du  premier  ordre 


et,  par  suite,  en  la  différentiant,  celle  du  second  ordre 


l’on  a,  non  seulement  cp ( Ç4  ) = 1 , mais  aussi  cp"(^  ) = — : et  les 

deux  formules  (44)  deviennent 


* X, 

2 ac  rs 


a — oj  — 


3 ji./nac  X. 

16  y 2 a ’ 


d’où 


m ( a — tu  ) — 


3 TC  IA2  X2 
1 28  a2 


(5i) 
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En  remplaçant,  dans  les  deux  premières,  c par  le  troisième 
membre  de  (47)5  puis  comparant  les  résultats  aux  expressions 
de  m et  de  a — to  qui  résultent  de  (38),  on  voit  que,  pour  l’intu- 

isolée,  le  coefficient  d’extinction  m égale  la  fraction  — > 


mescence 


3 y At 


ou  environ  les  0,886,  et,  le  retard  a — w,  les  -X— = o,665,  de 


leurs  valeurs  respectives  dans  un  système  d’ondes  périodiques 
dont  chacune  aurait  sa  moitié  convexe  de  même  hauteur  et  de 
même  énergie  que  l’intumescence  isolée  tout  entière. 

Quant  au  produit  m(a  — w),  indépendant  de  c comme  il  l’est 
de  n dans  le  cas  d’ondes  périodiques,  la  troisième  formule  (5i) 
donne  pour  sa  valeur,  à l’intérieur  d’une  conduite  circulaire  d’un 
diamètre  connu  D (d’où  <7  = ^ttD2),  et  en  se  rappelant 

finalement  que  p.  = o,oo58, 


(52) 


m ( a — o>  ) 


371^2 

8 D2 


. . . o,oooo3q6 

(environ)- — ■ — 2 ou 


0,00004 
D2  ’ 


Cette  formule  doit  être  approximativement  applicable  à toute 
intumescence  isolée;  car  une  telle  intumescence  aura  toujours  sa 
partie  antérieure  (de  la  tête  ou  du  front  jusqu’au  sommet),  qui  seule 
intervient  dans  les  relations  (44)^  plus  ou  moins  ressemblante  à 
celle  de  l’onde  dont  l’équation  est  (43). 

Avec  la  valeur  ni  — 0,0000866  obtenue  par  MM.  Violle  et  Vau- 
thier  pour  un  tuyau  de  om,7  de  diamètre  (où  ils  ont  trouvé  en 
outre  33om,  33,  après  réduction  à la  température  de  o°  et  à un 
état  hygrométrique  nul),  cette  relation  donne 

a — w = om,94, 

valeur  du  même  ordre  que  les  différences  respectives  des  vitesses 
de  propagation  constatées  par  ces  physiciens  ou  antérieurement 
par  Régnault  dans  des  tuyaux  de  divers  calibres. 

On  en  déduit  a = 33om,  33  -h  om,  g4  = 33 1 m,  3 environ,  pour  la 
vitesse  du  son  dans  l’air  libre,  sec  et  à la  température  de  o°C. 


XVI.  Si  l’intumescence  que  représente  à l’époque  t0  l’équa- 
tion (45)  conservait  indéfiniment  cette  équation  à part  l’introduc- 
tion, à son  second  membre,  du  facteur  lentement  décroissant  qui 
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exprimerait  sa  hauteur  h , son  énergie  totale  proportionnelle  à 
f ^ dx,  varierait  comme  h 2 ou  comme  une  exponentielle  de  la 
forme  . Le  coefficient  d’extinction,  que  j’appellerai  m\ 

propre  à cette  énergie  C,  et  dont  l’expression  à l’époque  t est  évi- 
demment 


(53) 


1 — ~ 1 cl_  r 2 7 
C a dt  a j de  dl  ^ C X ' 


aurait  donc  la  valeur  am  (*).  Mais  la  condensation  y décroît  plus 
vite  au  sommet  qu’elle  ne  fait,  du  moins  en  général,  aux  autres 
régions  de  l’onde,  ou  autrement  dit,  l’intumescence  s’allonge  tout 
en  devenant  moins  aiguë,  de  manière  à conserver  relativement  plus 
d’énergie  que  de  hauteur  ; et  le  rapport  de  m'  à ‘j.  m est  notablement 
inférieur  à l’unité. 

Divisons,  en  effet,  par  la  première  (5i),  la  moitié  du  der- 
nier membre  de  (53),  où  nous  savons,  d’après  (46  bis),  que 


(*)  C’est  justement  ce  qui  a lieu  clans  le  cas  des  mouvements  pendulaires  où 
la  formule  de  cp  est  (36).  Alors  le  dernier  terme  de  l’équation  (35),  qu’on  peut 


se  simplifie,  à cause  de  la  valeur  de  cp'(^+  (3S),  savoir 


inn  /2  met 

— sin  

a \ a 


2 mz  f . 2/nrS 

sin cos 

a \ a 


2 711TÇ 


— - — h cos 
a a 


qui  donne 


P 00  ,.y  70  tz  / n(  . 2/ixl- 

Jt , + — +cos  — 

Le  produit  de  cette  dernière  expression  par  cp(l*), 

iz  / n ( ï . 4 tit-X 

i / - - sin  -f-  cos2 > 

2 y a\2  a a ) 

peut  évidemment,  sous  le  signe  d’intégration  relatif  à être  remplacé  par  sa  va- 
leur moyenne;  ce  qui  permet  d’en  supprimer  le  premier  terme  entre  parenthèses, 

à valeur  moyenne  nulle,  et  il  reste  — Ainsi,  le  dernier  terme  de 

l’équation  (35),  valeur  de  la  dérivée  en  t de  j^dx,  devient 

— \x.\Jmz  |/«p(Ç)*rfÇï 

et  le  dernier  membre  de  (53),  expression  du  coefficient  m',  est  lui-même 
^ ^-nJl  h , justement  le  double  de  la  valeur  (38)  de  m. 


a 
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/tp(£)2a?£  = fitc  et  où  la  dérivée  en  t de  f y - dx  est  réductible  au 
dernier  terme  de  (35).  En  tenant  compte  de  la  troisième  rela- 
tion (46)  et  de  la  formule  (49)  ou  plutôt  de  celle  qu’en  a déduite 
une  différentiation  en  £,  il  viendra,  toutes  réductions  faites, 


Or  les  deux  intégrales  eulériennes  figurant  au  troisième  membre 
s’expriment,  par  la  formule 


B (p,  q ) = 


T(p)T(g) 

r {p  + q)' 


au  moyen  de  valeurs  de  la  fonction  T que  permettent  ensuite  de 
réduire  les  relations,  bien  connues, 

r(4)=i.2.3,  T(n)  = (n  — i)T(n  — 1),  r(/i)  T(i  — n)  = — — 


Ces  deux  intégrales  deviennent  ainsi,  respectivement,  — -)-?  --  , 

64  y 2 64  y 2 

et  la  formule  (54)  donne  enfin 

(55)  - = ——  =0,1  environ,  ou  m!  — m\/i  (1). 

2m  ^ 


(’)  Des  calculs  trop  longs  pour  trouver  place  ici  m’ont  permis  de  reconnaître 
que  le  coefficient  d’extinction  analogue,  m" , propre  à l’énergie  de  la  partie 
antérieure  de  l’onde  (depuis  le  front  jusqu’au  sommet),  est 


(56) 


:V/Î 


log(i  -f-  v/2)j  m = (v^  4-  °,374)/rc  ou  1.788m. 


La  partie  antérieure  s’allonge  donc,  mais  bien  moins  que  l’ensemble.  Si  l’on 
admet  d'ailleurs  que  son  équation  se  maintienne  assez  analogue  à (45)  pour 
permettre  d’y  regarder  à toute  époque  l’intégrale  f^2dx,  d’après  (46  bis),  comme 
le  produit  du  carré  de  sa  hauteur  h par  -‘tu c,  les  formules  (5i)  et  (56)  conduiront 
assez  facilement  aux  lois  de  variation  de  c et  h exprimées  par  la  double  égalité 


(67) 


Tcïc  - \jiac^  4-0,106 


ulv/tu  7 , 

— zat  = \2aco 

4 


c0  et  h0  étant  les  valeurs,  censées  données,  de  c et  de  h à l’époque  t — o. 
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Donc  l’onde  décroît,  proportionnellement,  moins  en  énergie 
qu’en  hauteur  (condensation  maxima);  ce  qui  exige  un  accroisse- 
ment de  sa  longueur  sensible  propre  à y compenser  en  partie 
V abaissement  du  maximum.  On  doit,  en  effet,  pouvoir  attribuer 
approximativement  à l’intumescence  une  longueur  sensible  totale 
(tant  à l’arrière  qu’à  l’avant  du  sommet),  proportionnelle  au  quo- 
tient de  l’énergie  par  le  carré  du  maximum  de  y.  Or  cette  longueur 
totale , ainsi  appréciée,  variera  ici  comme  rexponentielle 

g(2  m-rn')a(t—t0)=  e^2—\/V  ma(t—t0)  — g0,586  /na{t—t0)m 

Par  conséquent,  l’onde  paraît  s’allonger,  mais  moins  vite  qu’elle 
ne  s’abaisse  ( 1 ). 


(*)  La  réussite,  quand  cp(i-)  a l’expression  (45),  des  intégrations  indiquées  au 
dernier  terme  de  (35),  conduit  à déterminer  de  même  approximativement  l’ex- 
tinction d’une  onde  solitaire  par  la  voie  suivie  dans  mon  Mémoire  cité  de  1878 
( Journal  de  Mathématiques,  t.  IV,  p.  36a  à 364),  mais  en  attribuant  à cette 
onde  l’équation  même  (45),  au  lieu  de  l’équation  exacte,  presque  analogue, 

t / ] j 

cp  (£)  = coh-2  -,  avec  c = 2 H 4/  ^7  et  a = t*>  = \J g H. 


Alors  on  introduit  aisément  dans  l’équation  (42)  du  Mémoire  nos  variables 
actuelles  £ et  rt,  à la  place  de  t et  n,  en  observant  que  £ = x — at,  p = n\[â , 
et/(0  = ?(£)• 

Cette  équation  (42)  devient,  grâce  à l’eflectuation  finale  des  calculs  comme 
ci-dessus  et  après  élimination  de  a , c, 


On  en  déduit,  pour  le  coefficient  e'  qui  figure  aux  équations  de  la  p.  364, 
valeur  ^ \3tzs.  \/  p2  \Jg,  et  la  formule  définitive  (47)  donnant  l’énergie  décroissante 

O 

C de  l’onde  est 


1 

76 


dx 


9 i 

H U* 


(Extrait  du  Journal  de  Physique , 2e  série,  t.  X;  juillet  1891.) 
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